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© A two channel real time octave analyzer is equipped with an adaptive time domain phase compensation filter 
whose poles and zeros are programmably selected to counteract phase error associated with any given pair of 
microphones. Precise matching of microphone phase characteristics is thereby achieved. In the preferred 
embodiment, the transfer function for the adaptive filter is determined using a pole/zero synthesis technique 
based on cross spectrum data acquired in an FFT spectral analysis of the probe microphones. 
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Field of the Invention 



The present invention relates to a method and apparatus for precisely matching the characteristics of a 
pair of transducers, such as microphones, for measurement applications. 

5 

Background and Summary of the Invention 

For expository convenience, the present invention is illustrated with reference to the matching of 
characteristics of microphones used in a probe for measuring sound intensity. It should be recognized, 

io however, that the invention can advantageously be applied in a number of other situations requiring precise 
matching of phase and/or amplitude characteristics. 

Sound intensity is a vector measurement of the average rate of sound energy transmitted in a specified 
direction through a unit area normal to this direction at a specific point. Such measurements are commonly 
used to quantify, for example, the noise emitted from industrial equipment or machinery. 

75 Sound intensity measurements are conventionally made using a pair of precisely matched, closely 

spaced microphones (commonly known as a sound intensity "probe"). As explained more fully in Pope, J., 
"The Two-Microphone Sound Intensity Probe," Journal of Vibration, Stress, and Reliability in Design, Vol. 
110, January, 1988, pp. 97-103, sound intensity is related to the cross-spectrum detected by the two 
microphones in the frequency domain according to the following formula: 



\ r (uy) = -lm(G AB ypwAr (1) 

The equivalent expression in the time domain is: 



I r = ~(Pa + Pb)/(Pb * P.)<*t (2) 

2pAr 

30 . . 

where: 

p is the density of the acoustic medium; 
oo is the frequency in radians; 
Ar is the effective microphone separation; 
35 p a is the sound pressure at microphone a; 
p b is the sound pressure at microphone b; 
l r is sound intensity in direction f from a to b; 
and 

lm(G A B) represents the imaginary part of the cross-spectrum between p a and p b - 
aq The measurements necessary to make the , above-detailed calculations can be performed using two 
different types of instruments: Fast Fourier Transform (FFT) spectrum analyzers and real time octave 
analyzers. Each has its respective advantages and disadvantages. 

An FFT analyzer, such as the Hewlett-Packard HP 35660A and 35665A, operates by digitally sampling 
an analog input signal and performing a fast Fourier transform on the sampled data to determine its spectral 
45 composition. The results of the Fourier analysis are a series of spectral coefficients, one corresponding to 
each of a plurality of frequency "bins." 

A drawback of the FFT analyzer technique is that the frequency bins into which the Fourier analysis 
resolves the spectral composition are uniformly spaced in frequency. A 400 bin analysis of the spectrum 
between 0 and 10 KHz, for example, results in bins that each correspond to 25 Hz. While such resolution is 
so more than adequate for higher frequencies, it is inadequate at low frequencies. 

t To provide adequate resolution at low frequencies, a second FFT measurement spanning a smaller 
frequency range (such as 0 - 1 KHz) is generally required. The results of the two measurements are then 
combined to yield a final result. This procedure, however, is problematical since a single instrument cannot 
make both measurements at the same time. The combined measurement instead reflects two different 
55 measurements made at two different times. This non-real time operation is often unacceptable. 

Real time octave analyzers, in contrast, use a series of bandpass filters (often implemented digitally in a 
sampled data system) to determine spectral composition. These filters are generally centered at 
logarithmically-spaced frequencies (often in one-third octave steps), thereby providing increasingly finer 
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resolution at increasingly lower frequencies. 

(The human ear perceives sound in a logarithmic fashion, making octave-based analysis a popular 
mode of measurement. Further, most acoustic standards, such as the IEC 1043 Instruments for the 
Measurement of Sound Standard, are specified in this fashion.) 
5 Several real time octave analyzers are known, including the Bruel & Kjaer 2133 and the HP 35665A. 

(The HP 35665A, which was noted above as an FFT analyzer, combines FFT and real time octave analysis 
capability in a single instrument.) In the B&K 2133, sound intensity is calculated by combining the outputs 
of the bandpass filters to perform a sampled data counterpart to the operation described in equation (2) 
above. 

io In order for the cross-channel data acquired in either an FFT or real time octave analyzer measurement 

setup to be useful, the characteristics of the probe microphones must be precisely matched. Amplitude 
matching is relatively easy to obtain. (Further, the sound intensity calculation is less sensitive to amplitude 
mismatch. Cf. pope, supra.) Phase matching, however, is difficult. 

Two primary mechanisms contribute to phase mismatches between otherwise identical microphones. 

75 The first relates to variability in diaphragm damping. This phenomenon is particularly evidenced at 
frequencies above a few hundred Hz. Due to manufacturing tolerances, otherwise identical microphones 
may commonly exhibit a difference in phase response of several degrees when measured at one kilohertz. 

The second 1 mechanism that contributes to microphone phase mismatch is an interaction between a 
cavity behind the microphone diaphragm and a hole venting this cavity to the atmosphere, effectively 

20 creating an acoustical low pass filter. Differences in this static pressure equalization between microphones 
can cause a phase variance between otherwise matched microphones of up to three degrees at twenty 
hertz. (Twenty hertz is typically the low frequency limit for sound intensity measurements.) 

In order for a pair of microphones to be successfully utilized for sound intensity measurements, their 
phase characteristics should generally be matched to within less than 0.3 degrees throughout the 

25 measurement span of interest. (This 0.3 degree figure is highly application dependent. In some instances a 
variance of up to 1 degree may be acceptable. In others, a variance of less than 0.1 degrees may be 
required.). As noted, current manufacturing processes can only reliably match microphones to within about 
3 degrees, an order of magnitude higher than this approximate threshold requirement. Accordingly, 
microphones must generally be matched by a tedious manual selection process. 

30 Low frequency phase error is usually the most critical in making sound intensity measurements. As 
noted, sound intensity is a vector measurement. The direction of the measured sound is determined by the 
phase delay with which the same sound reaches the two microphones. Since, at low frequencies, the sound 
wavelength is several meters long, and the microphones may only be spaced by a centimeter or two, the 
cos phase delay associated with off-axis sound vectors is quite small and would be masked by even 

35 small phase errors (<f> is here the angle between the intensity vector and the center line of the microphones) 

If the problem of low frequency phase error could be solved, the remaining task of selecting 
microphones by manually matching high frequency responses would be considerably simplified. 

While in real time octave analyzers, the microphones themselves must be precisely matched, in FFT 

40 analyzers it is possible to match the microphones using correction factors internal to the analyzer. In 
particular, it is possible to determine the phase error between a pair of microphones in the frequency 
domain, and then multiply the cross spectrum between the microphones by the complex conjugate of this 
error term to effect correction. This method has been used for many years in FFT instruments. 

(To determine the phase error between a pair of microphones, a number of different methods may be 

45 employed. Exemplary are those disclosed in Chung, J.Y., "Cross-Spectral Method of Measuring Acoustic 
Intensity Without Error Caused by Instrument Phase Mismatch," Journal of the Acoustical Society of 
America, Vol. 64, No. 6, 1978, pp. 1613-16, and Seybert, A.F., "Measurement of Phase Mismatch Between 
Two Microphones," NOISE-CON 85 Proceedings, 1985, pp. 423-28.) 

Real time octave analyzers have not previously been susceptible to this type of phase correction. 

so Consequently, extensive research has been conducted in methodologies of producing closely matched 
microphones. Exemplary of this work are: Frederiksen, E., "Phase Characteristics of Microphones for 
Intensity Probes," Proceedings of 2nd International Congress on Acoustic Intensity, Senlis, 1985, pp. 50-57; 
and Frederiksen, E., et al, "Pressure Microphones for Intensity Measurements with Significantly Improved 
Phase Properties," Bruel & Kjaer Technical Review, No. 4, 1986, pp. 11-21. The latter paper proposed a 

55 solution to low frequency phase mismatch by mechanical compensation -- namely the addition of two 
additional cavities behind the diaphragm to attenuate the effect of low frequency vent-cavity resonances, 
and hence minimize phase variabilities between microphones. Microphones optimized for phase accuracy 
are also shown in U.S. Patents 4,887,300 and 4,777,650. While effective, such approaches result in 
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significantly increased manufacturing costs. 

In accordance with the present invention, a real time octave analyzer is equipped with an adaptive time 
domain phase compensation filter whose poles and zeros are selected to counteract the low frequency 
phase error associated with any given pair of microphones. Precise equalization of unmatched microphone 
5 phase errors is thereby achieved without resort to elaborate mechanical compensation schemes. In the 
preferred embodiment, the transfer function for the adaptive filter is determined using a pole/zero curve 
fitting technique based on cross spectrum data acquired in an FFT spectral analysis of the probe 
microphones. 

The foregoing and additional features and advantages of the present invention will be more readily 
w apparent from the following detailed description, which proceeds with reference to the accompanying 
drawings. 

Brief Description of the Drawings 

75 Fi 9- 1 is a schematic block diagram of a real time one-third octave analyzer according to one 

embodiment of the present invention. 

Fig. 2 is a plot showing a cross spectrum of pink noise sampled by a pair of unmatched microphones in 
a reference state and in a switched state as a function of frequency. 

Fig. 3 is a plot of the microphone phase error determined from Fig. 2, and a correction curve 
20 synthesized with three poles and three zeroes. 

Fig. 4 is a plot of the microphone phase error after correction by the apparatus of Fig. 1. 
Fig. 5 is a block diagram of a canonic direct form II implementation of a time domain phase correction 
filter used in the analyzer of Fig. 1 . 

25 Detailed Description 

Referring to Fig. 1, a one-third octave analyzer 10 according to the present invention includes two 
measurement channels A and B. Only channel A is illustrated. Channel B is identical except for the 
omission of a filter 12 detailed below. 
30 Channel A includes a microphone 14, an analog anti-alias filter 16, an analog-to-digital (A/D) converter 
18, and a series of analysis stages 20. Each analysis stage 20 includes a digital anti-alias low pass filter 22, 
a decimation stage 24, and three one-third octave digital bandpass filters 26/26726". 

The microphone 14 is one of two that comprise the sound intensity probe. The output from this 
microphone is low-pass filtered by the analog anti-alias filter 16 to attenuate spectral components above 
35 25.6 KHz. This band-limited signal is then sampled 65,536 (2 16 ) times per second by the A/D converter 18. 

The sampled signal from the converter 18 is provided to the first analysis stage 20a, which low pass 
filters the incoming sampled data with filter 22a to attenuate components above 12.8 KHz. The filtered 
signal is then decimated by decimator 24a. (Decimation refers to a known process wherein periodic data 
samples not essential to fulfillment of the Nyquist criterion are ignored. Here, decimation permits all lower 
40 full-octave bands to be analyzed with the computational power that would normally be required to analyze a 
single octave at the original 65.536 KHz sample rate.) 

The decimated output signal (now consisting of 32,768 samples/second) is applied to the three digital 
bandpass filters 26a, 26a 1 and 26a" (ANSI Std. S1.11 1986 order 3, type 1-D filters), which resolve the 
sampled spectra into one-third octave bands centered at 8000, 6250, and 5000 Hz, respectively. 
45 The output from the first decimator 24a provides the sampled data used by the second analysis stage, 
20b. This data is again low pass filtered (this time at 6.4 KHz), decimated (to 16,384 samples/second), and 
provided to one-third octave digital filters 26b, 26b' and 26b" centered at 4,000, 3,125, and 2,500 Hz 
respectively. 

Fortuitously, since decimation halved the sample rate, and since the analysis bands are also related by 
so two (per ANSI Std. S1. 11-1986, base 2 system of digital one-third octave filters), the positions of the poles 
and zeros of filters 26b/26b726b" are the same as those of filters 26a/26a726a". This greatly facilitates 
implementation of the bandpass filters since the same set of three filters used in stage 20a is simply 
replicated in stage 20b. - . 

Additional analysis stages 20c, 20d and 20e are similarly cascaded, one after another. Each low pass 
55 filters the input signal at half the prior filter frequency, decimates by two, and resolves the decimated data 
into the next three lower one-third octave bands. Again, decimation of the data, coupled with the halving of 
the analysis bands, permits the same set of three filters to be used in each successive stage. 

Analysis stage 20f includes a filter 12 not found in the other analysis stages, nor found in the other 
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measurement channel ("B"). The illustrated filter 12 is an adaptive, unity gain, time domain, infinite impulse 
response (MR) filter whose purpose is to counteract phase errors associated with the probe microphones 14 
at frequencies below 300 Hz. The particular manner in which the poles and zeros of this filter are 
determined is detailed below and illustrated with reference to Figs. 2-4/ 
s First, the phase error between the two probe microphones is determined. A suitable method is 

described by Chung, supra. Chung's method basically is a circuit switching technique wherein the cross 
spectrum is measured once (yielding a measured spectrum Gab) and then measured again with the 
microphones switched (yielding a measured spectrum G s A b)* The switching permits the phase response due 
to actual intensity to cancel upon division, leaving just the probe phase error term in the following equation: 

10 

e j " e = (G u / G s \b)* (3) 



75 where 



20 is the phase error as a function of frequency; 
and 



is the complex conjugate of G S AB - - '■ • - 

Similarly, the phase error cancels upon the multiplication of Gab and 

30 ' s * 

AB 

yielding an estimate of the cross-spectrum which is devoid of errors due to the phase mismatch. This 

estimate would then be inserted in equation (1). 
35 Fig. 3 shows, by the curve 30, the conjugate result of this calculation for the G A b and G s A b data shown 

in Fig. 2. The conjugate of the error phase is used as a basis for the correction filter 12. 

It is possible to synthesize a correction filter transfer function that is arbitrarily close to curve 30. Curve 

32 in Fig. 3 shows a transfer function synthesized using three poles and three zeros. Such an approximation 

is adequate for most applications. 
40 The methodology used to synthesize a time domain filter transfer function from data such as that 

represented by curve 30 has been known in the art for many years and is illustrated, for example, in U.S. 

Patents 4,654,808, 4,654,809 and 4,658,367, the disclosures of which are incorporated herein by reference. 

In the best mode, the Chung microphone phase error measurements are made on a HP 3563A FFT 

analyzer that has the built-in capability to fit a synthesized S-domain transfer function to an arbitrary curve 
45 using the methodology disclosed in these patents. 

The output of the curve fitting process are data representing the poles and zeros of the synthesized 

transfer function. 

In some instances, the curve-fitting methodology described in the above-mentioned patents and 
implemented in the HP 3563A analyzer may result in poles at the right hand side of the S-plane, making the 
50 resulting filter unrealizable. To avoid this result, it has been found helpful to fix, a priori, one of the poles at 
0 Hz (DC). (Further a priori placement of poles and zeros can aid in generating uniform gain factors from 
the curve fitting process. A physical justification for choosing additional fixed pole locations may be found in 
the cutoff frequency of the acoustical low pass filter created by the microphone vent-cavity interaction.) 

In the illustrated embodiment, with one of the poles constrained to the origin of the S-plane, the 
55 pole/zero synthesis function yielded the following poles and zeros (each of which here lies lie along the real 
S-axis): ' 
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POLES 


ZEROS 


-135 


-0.8962 

-6.0 
-132.33 



The HP 3563A normally determines the synthesized transfer function in the S plane. However, 
conversion into the Z-domain is straightforward using the bi-linear transform - a feature which is also 
implemented as an automated operation in the HP 3563A analyzer. 

Using a pre-warp frequency of 200 Hz, the following pole-zero locations are obtained in the Z-domain: 



75 



POLES 


ZEROS 


1.0 
.355575 
.992283 


.958618 
.364269 
.993708 



(Pre-warp refers to the frequency at which the S- and Z-domain representations coincide.) 
20 The filter transfer function in the Z-domain is represented by: 



Afftl-C/ 1 ) f>,z~ r 

r- 1 r-0 

25 H(z)= —. = (4) 

3 ' 3 

IKl-dfcZ' 1 ) 1- Ia k z' k 

k-l' k-1 

30 . . 

The pole-zero equation for H(z) can be written by setting c, equal to the zero locations and d k equal to 
the pole locations. Once in the Z-domain, the transfer function can be expanded into polynomial form. The 
polynomial coefficients of z can then be used to implement the synthesized correction filter in a 
configuration called canonic direct form II. Non-unity gain of the correction filter is scaled by a factor of A. 

35 Such an implementation for the above-described filter is shown in Fig. 5 and will be recognized as an 
infinite impulse response (IIR) topology. In this example, the coefficients are real, although they need not 
be. The coefficients are also desirably reprogrammable, permitting the filter to be adapted to different 
microphones, or to track errors in a single set of microphones over time. 

Fig. 4 illustrates the probe phase error after correction by the time domain correction filter 12 (curve 34) 

40 and compares it with the phase error prior to correction (curve 30). As can be seen, the phase errors are 
greatly minimized. , 

(The curves of Figs. 2-4 were obtained in a proof of concept measurement using randomly selected 
microphones in a non-anechoic environment. Error correction substantially in excess of that portrayed can 
be obtained in more controlled conditions.) 

45 Returning to Fig. 1, the output from the correction filter 12 is provided to one-third octave filters 
26f/26f726f\ which are centered at 250, 200 and 160 Hz, respectively. A cascaded arrangement like that 
described above then follows, with analysis stages 20g and 20h producing outputs for the bottom-most two 
octaves (125/100/80 Hz and 62.5/50/40 Hz). 

To determine sound intensity, the outputs from the one-third octave bandpass filters 26/26726" of 

so channel A are combined with their counterparts of channel B in the conventional manner expressed by 
equation (2). 

In the preferred embodiment, the analysis stages 20, including the time domain correction filter 12, are 
implemented using a general purpose digital signal processing circuit, such as the Motorola MC56001. In an 
alternative embodiment, the bandpass filters alone are implemented using a general purpose DSP circuit 
55 and the decimators 24 and low pass filters 22 of the analysis stages are implemented in a custom gate 
array. 

In applicants best mode, the above-described steps of microphone calibration and correction are 
integrated in an automated procedure that is programmed into a multi-channel FFT/real time octave 
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analyzer instrument. Upon execution of this procedure, the instrument: provides microphone setup instruc- 
tions to a user, provides a source of pink noise, acquires fin FFT mode) a first cross-spectrum, instructs the 
user to switch the microphones (assuming the Chung procedure is employed), acquires a second cross- 
spectrum (again in FFT mode), performs the operation of equation (3) to determine the phase error, 

5 synthesizes coefficients for a correction filter in accordance with the patented curve-fitting procedure, 
implements the correction filter in the real time octave analyzer architecture of Fig. 1 , and then indicates to 
the user that the system has been calibrated and compensated and is ready to make sound intensity 
measurements. Data that is thereafter acquired by the calibrated system is processed in accordance with 
the sampled data counterpart of equation (2) to yield a real time sound intensity measurement. 

w Having described and illustrated the principles of my invention with reference to a preferred embodi- 
ment thereof, it will be apparent that the invention can be modified in arrangement and detail without 
departing from such principles. For example, while the invention has been illustrated with reference to a 
time domain correction filter optimized just for low frequencies, it will be recognized that the same 
principles permit compensation over an arbitrarily broad bandwidth. (More than three poles and zeros may 

75 be desirable in filters adapted to compensate the entire measurement band.) Similarly, while the invention 
has been illustrated with reference to an apparatus employing a single adaptive time domain correction 
filter, it will be recognized that in other embodiments, a plurality of such filters may advantageously be 
employed. Multiple serial correction filters may permit more precise tailoring of the phase correction curve. 
One filter may correct for higher frequency phase error (and have nil effect at low frequencies), and another 

20 filter may correct for lower frequency phase error. Alternatively, multiple filters may be used by interposing 
one between the decimator 24 and the filters 26/26726 ,, of each analysis stage 20, rather than serially in the 
chain of analysis stages. Still further, multiple filters may be used by placing one or more in each 
measurement channel, in which case it is the phase difference between the filter transfer functions that is 
effectively implemented. 

25 While the invention has been illustrated with reference to a three pole/zero correction filter, it will be 
recognized that other filters may be found appropriate in other situations. For example, if the microphones 
are fairly closely matched, then a single term correction filter may suffice. Similarly, while the invention has 
been illustrated with reference to an S-domain curve fitter with a subsequent bi-linear transform to generate 
the 2-domain topology, in other embodiments the poles and zeros can be determined directly by a Z- 

30 domain curve fitter, such as that found in the HP 3563A. Finally, while the invention has been illustrated with 
reference to an adaptive MR filter topology for the time domain correction filter, it will be recognized that 
other phase compensating filter topologies may readily be employed. 

In view of the wide variety of embodiments to which the principles of my invention can be applied, it 
should be apparent that the detailed embodiment is illustrative only and should not be taken as limiting the 

35 scope of my invention. Rather, I claim as my invention all such modifications as may come within the scope 
and spirit of the following claims and equivalents thereto: 

Claims 



40 1. In a method of time domain spectral analysis, an improvement comprising the steps: 
determining a transfer characteristic sought to be compensated; 

synthesizing a pole/zero transfer function approximating said transfer characteristic; and 
programming an adaptive digital filter in accordance with said synthesized transfer function. 

45 2. A method of determining sound intensity comprising the steps: 

determining a phase error function associated with signals from first and second microphones; 
processing ah output signal from the first microphone through a first series of bandpass filters to 
determine the spectral composition thereof; 

processing an output signal from the second microphone through a second series of bandpass 
so filter to determine the spectral composition thereof; and 

analyzing output signals from said first and second series of bandpass filters to determine sound 
intensity; 

wherein at least one of said processing steps further includes time domain filtering to correct for 
said determined phase error function. 



55 



The method of claim 2 which further includes: 

synthesizing a transfer function approximating the conjugate of the determined error function; and 
programming a digital filter to implement said synthesized transfer function. 
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4. The method of claim 3 which further includes constraining a pole of the synthesized transfer function to 
DC when synthesizing said transfer function. 

5. The method of claim 2 which further includes using a single measurement instrument to: 
5 perform FFT analyses of signals from the first and second microphones; 

determine from said FFT analyses a phase error function associated with said microphone signals; 
synthesize a transfer function approximating the determined error function; 

adapt a programmable time domain digital filter to implement said synthesized transfer function; 

and 

70 compensate a digital data stream associated with the first microphone using said time domain 

digital filter. 

6. The method of claim 5 which further includes using said same measurement instrument to perform FFT 
analyses of cross spectra using the microphones in reference and switched states, and determining 

75 from said cross spectra a phase error function associated with the microphone signals. 

7. A real time octave analyzer having at least two measurement channels, each measurement channel 
comprising: 

a source of input signals; 
20 an anti-alias filter having an input coupled to the source of input signals; 

an analog-to-digital converter having an input coupled to an output of the anti-alias filter; 
a plurality of successive analysis stages, each having an input, an output, and at least one 
bandpass filter associated therewith; 

the input of a first analysis stage being coupled to an output of the analog-to-digital filter; 
25 the input of each successive analysis stage being coupled to an output of a previous analysis 

stage; , 

one of said measurement channels further including a time domain digital filter adapted to 
compensate for a phase error associated with the source of input signals. 

30 8. The analyzer of claim 7 in which the time domain digital filter is an infinite impulse response filter. 

9. The analyzer of claim 7 in which the time domain digital filter has a pole at DC. 

10. The analyzer of claim 7 in which the time domain digital filter is programmable so that it may be 
35 adapted to compensate for phase errors associated with a variety of input signal sources. 
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@ Time domain compensation for transducer mismatch. 



© a two channel real time octave analyzer is 
equipped with an adaptive time domain phase com- 
pensation filter (12) whose poles and zeros are prog- 
rammably selected to counteract phase error asso- 
ciated with any given pair of microphones. Precise 
matching of microphone phase characteristics is 



thereby achieved. In the preferred embodiment, the 
transfer function for the adaptive filter is determined 
using a pole/zero synthesis technique based on 
cross spectrum data acquired in an FFT spectral 
analysis of the probe microphones. 
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